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ABSTRACT 

The  energy  method  is  applied  to  study  the  stability  of  two 
types  of  difference  approximations  to  parabolic  partial  differ- 
ential equations,  the  alternating  direction  methods  of  Douglas, 
Peaceman,  and  Rachford,  and  a  new  semi-explicit  method.   Each 
difference  scheme  is  proved  to  be  unconditionally  stable.   These 
results  apply  to  parabolic  equations  with  variable  coefficients, 
defined  in  cylindrical  domains  with  an  essentially  arbitrary 
bounded  base. 
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ALTERNATING  DIRECTION  AND  SEMI- EXPLICIT  DIFFERENCE 
METHODS  FOR  PARABOLIC  PARTIAL  DIFFERENTIAL  EQUATIONS 


1 .   Introduction. 

In  previous  papers  [12],  [l'j>] ,  [l4]  the  author  developed  a 
difference  analogue  of  the  energy  method  for  determining  the 
stability  of  difference  approximations  to  partial  differential 
equations  with  variable  coefficients.   The  purpose  of  this  paper 
is  to  apply  this  method  to  establish  the  unconditional  stability 
of  two  types  of  difference  approximations  to  parabolic  differential 
equations,  the  (implicit)  alternating  direction  methods  of  Douglas, 
Peaceman,  and  Rachford  [;>],  [5],  [15] »  and  a  new  semi-explicit  method, 

For  the  model  problem,  the  first  boundary  value  problem  for  the 
heat  conduction  equation  in  a  rectangular  domain,  the  unconditional 
stability  of  the  alternating  direction  methods  was  proved  in  [3] 
and  [5]»   The  proof  consists  in  showing,  with  the  aid  of  Fourier 
analysis,  that  the  von  Neumann  stability  condition  [4],  [11]   is 
always  satisfied.   It  can  be  shown  [1],  however,  that  this  method 
of  proof  cannot  be  extended  beyond  the  model  problem. 

With  the  aid  of  the  energy  method  we  prove  that  the  results 
in  [~j>]    and  [5]  can  be  extended  beyond  the  model  problem^   We  first 
treat,  as  a  typical  case,  the  heat  conduction  equation  in  a  cylindri- 
cal domain  with  an  essentially  arbitrary,  bounded  base.   Then  we 
indicate  briefly  the  extension  to  parabolic  equations  with  variable 
coefficients . 

The  second  type  of  difference  method  we  term  the  semi-explicit 
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method,  because  it  is  an  explicit  method  only  for  certain  orderings 
of  the  net  points.   The  idea  for  this  difference  method  comes  from 
the  observation  that  there  is  a  formal  correspondence  between  para- 
bolic difference  equations  and  iterative  methods  for  solving  elliptic 
difference  equations;  the  semi-explicit  method  corresponds  to  the 
well  known  method  of  successive  displacements  [?]• 

The  only  other  known  example  of  an  unconditionally  stable 
explicit  difference  method  is  due  to  DuPort  and  Frankel  [6].   Their 
method,  however,  requires  two  lines  of  initial  data  to  start  the 
solution,  while  the  semi-explicit  method  is  self-starting.   On  the 
other  hand,  both  of  these  methods  involve  a  similar  local  truncation 
error,  and  they  must  be  subjected  to  a  mild  mesh  ratio  condition  in 
order  to  be  consistent  [11]  with  the  differential  equation  being 
approximated. 

For  other  applications  of  the  energy  method  to  the  stability 
problem  for  partial  difference  equations,  see  Friedrichs  [8], 
Kreiss  [9],  Lax  [10]  and  Lees  [12],  [13],  [l2*]. 

In  a  paper  to  follow,  the  energy  method  will  be  applied  to  the 
problem  of  determining  the  rate  of  convergence  of  iterative  methods 
for  solving  elliptic  difference  equations.   It  will  include,  in 
particular,  a  solution  to  this  problem  for  the  alternating  direction 
methods  relative  to  an  arbitrary  domain. 

2  .   Notation  and  Definitions. 

In  this  section  we  describe  certain  preliminary  concepts,  neces- 
sary for  our  formulation  of  the  difference  schemes.   We  denote  by 
ft  a  bounded,  open  subset  of   E„,   with  boundary   ft.   The  Euclidean 
length  of  a  point  x  =  (x1,  x2,...,  xN)  e  en  will  be  denoted  by 
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|x[.   Let  h  =  (h-,  ,  ho,...,  h„)  e  EN  have  positive  coordinates, 

and  define   G,   to  be  the  set  of  all  (net)  points  (iJi-^ighg, ..  ,iNhN)e  E] 

the  i  .   being  integers,  positive,  negative,  or  zero.   Two  points 

x>    y   £  G,   are  called  neighbors  if   |x  -  y|  =  h.,   for  some 

i,  i  =  1,  2,...,  N.   The  points  x   £    Gh  n   Q,     all  of  whose  neighbors 

belong  to  "Q,   the  closure  of  Q,,      we  denote  by   f^.   The  points 

x  e  Qh  -  £1   with  the  property  that  at  least  one  neighbor  belongs 

to  a       we  denote  by  i^.   Finally,  we  put  7^  =  %  u'  %• 

If   N  is  any  subset  of   EN,   we  define  c?(N)   to  be  the  collection 
of  all  real-valued  functions  defined  on   EN  whose  support   is  con- 
tained in   N.   Clearly,  ^(N)   is  a  real  linear  space  for  the  usual 
operations.   In  particular,  (^(fyj   is  a  real  finite-dimensional 
linear  space,  with  dimension  equal  to  the  number  of  points  in  Q,   . 

Let   S  =  [0,co]   and   S°  =  S  -  f 0 j  .   For  each  positive  number 
K,   the  time  step,  we  put 


SK 


=  [t  e  s|t  =  mK,  m  =0,  1,...,  ]■ 


and   S°  =  S„  H  S. 
i\  a. 

We  shall  approximate  the  solutions  of  differential  equations 
by  functions   t  — >  u(t),   defined  on   SK,  taking  their  values  in 

If   u  e  ^(ei\t)   we  define  the  linear  translation  operators   E— 
as  follows: 

The  support  of  a  function   f   is  the  closure  of  the  set 


[x|f (x)  t  oj 
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+1 
E   [u] (x)  =  u(x1, . . . ,  x±   +  h±, . . . ,    xN) . 

In  terras  of  these  we  define  the  first  order  forward  and  backward 
difference  operators: 

-lr^i 


and 


V±   u  =  h^[Ex  u  -  u] 


V  ±   u  =  h^tu  -  E  X   u] 


Evidently,  we  have  the  relation 


(2.1)        E-1^  u]  =V±   u, 


The  first  order  centered  difference  operator  V.      is  defined  by 


(2.2)        V.  u  =  |[V±  +  VJu  -  ^-[E1  -  E_1]u. 


Difference  operators  of  higher  order  are  defined  in  the  obvious 
way,  by  repeated  application  of  these  formulas. 

For  a  function   t  — >•  u(t)  e  £(EN)  >      defined  in   S   (resp.  S°) , 
we  use  the  notation 

u, (t)  =  K_1[u(t+K)  -  u(t)] 


and 


u  (t)  =  K  2[u(t)  -  u(t-K)]. 
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J>.      The   Linear   Space     £(fyj  - 

We  provide       §  (fy.)      with  an   inner  product   defined   as   follows 
if     u,    v   e    £(%), 


(u,    v)    =  hN       V"     u(x)v(x), 


x   e    Gh 


where  h   =  huhu...  h...   Associated  with  this  inner  product  is 


the  norm   ||u||  =  ^/(u,  u) 
The  formula 


N   _     2iV2 


i=i 


defines  another  norm  for  (?(%>)• 

Since  &({l)   is  finite-dimensional,  these  norms  are  equivalent; 
that  is,  there  exist  two  constants  m  and  m,   such  that 


H|u||2   <  ||u||2   <  mjulf, 


for  every  u  e  §(fy1)-   Of  course,   m  and  m.   depend  on  the  dimension 
of  ^(ft,)   and  therefore  on  the  (net  spacings)   h.. 

We  may  assume  that  TL       is  contained  in  the  rectangle 


R  :  a±  <  x±  <  b1, 


where   a.   and  b.   are  integral  multiples  of  h. 
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LEMMA  1.       If      u   €   SifyJ      then 


m|MI     <   lul^    , 


where  m  is_  the  minimal  eigenvalue  of  the  Laplace  difference 
operator 

N  _ 

(3.1)  Ah(u)  =y  Vi  V±  u 

i=l 


relative  to  the  net  region  ft,  ,      and 


"  >  4  XI  i       L2(b.-a.)J 


Proof.   The  minimal  eigenvalue  m  can  be  characterized  (Courant, 
Friedrichs,  and  Lewy  [2])  as  follows: 

l|u|g 

m   =  inf  -x   , 

0  fc  u   e  $(0^      ||u||2 


and   therefore 


m  < 


II  u  111 
l|u||2    ' 


for  any  u  e  <^(a),   not  identically  zero.   This  proves  the  first 
part  of  the  lemma. 

Now  let  R,  =  R°  O  Gu,   where  R°  is  the  interior  of  the 

h         h 
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rectangle     R.      Denoting  by     m(R.)      the  minimal   eigenvalue   of   A 
relative   to     R,  ,      we  have 


h 


/      x                                           1HI1 
m(R   )    =  inf  1 

h  Ofuc   <?(Rh)     Ijujp 


<  inf 


l|u|g 


0   fc  u   e   <J(fih)      ||u|r 


in. 


But   it   is   well   known    [7]    that 


N        _2  2f     irhj 


m(ig    -*|~  V   sin^^^-j. 


i=l 


and  this  completes  the  proof  of  the  lemma. 
LEMMA  2.   If   u  e  fc^)   then 


^ig  <  *<sr  hi2)  imi2- 


Proof.      We  have   that 
2.F}       ii2 


hJ|V.uf    =  h1Ny~   [u(x)    -    E"x[u](x)] 


x  e    G 


h 


<  2h 


N 


[u2(x)    +   E_i[u2](x)] 


x   e    g 


h 


=  *  iiuir 
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From  this   we   obtain 


l|u|g   =^[11^  u||2  <  4(^h-2)||u|p, 
i=l  i=l  -1 


which  is  what  we  set  out  to  prove. 

Remark :   From  lemma  1  we  see  that  m  is  bounded  away  from  zero, 

independently  of  h: 

m  >  4  J"  (b.  -  a±)  *. 

i=l 

On  the  other  hand,   m,   is  singular  at  h  =  0.   This  is  to  be 
expected  since  differentiation  is  an  unbounded  operator  in  Lp. 

4.   Difference  Methods. 

In  this  section  we  describe  the  difference  methods  for  the  first 
boundary  value  problem  for  the  parabolic  equation 

(4.i)  ^  =  y~  <UJ       (x  e  a,  t  >  o). 


The  function  t  —>   U(t)  e  &(ft) ,   for  t  >  0,   is  assumed  to  be 
sufficiently  smooth  and   U(t)  e  C  (0.)      for  a  suitable  integer   P, 

A  function  t  ->  u(t)  e  £(7^) ,   defined  for  t  e  SR  will  be 
called  admissible  if  u(0)  =  U(0),   and  u  =  U  on  J^  X  S°.   In 
general,   U  is  not  known  in  fi,  X  S°,   so  that  we  make  the 
assumption:  there  exists  a  null  sequence  lir    f     of  net  spacings 
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such  that   £18  c  h.      We  always  assume  that  h  belongs  to  the 
sequence   (h  (•   With  this  assumption,  an  admissible  function  is 
uniquely  specified  in  the  complement  of   11  X  S„   (relative  to  G,), 
in  terms  of  the  initial  and  boundary  values  of  U. 

In  the  semi-explicit  difference  method  for  (4.1),  one  determines 
an  admissible  function  t  — >  u(t)   such  that,   in  fyX  S„.,   the 
following  difference  equation  is  satisfied: 

N    ,  — 

(4.2)  ut  =Ah(u)  -  K^T"  h.1  V.  ut, 

i=l 

where  A,   is  the  Laplace  difference  operator  defined  in  (3«l). 
The  classical  explicit  difference  equation 

(4.3)  ut  =Ah(u) 

is  known  to  be  conditionally  stable  [7]*  the  condition  being  that 
the  mesh  ratio 


A  =K^"h.^ 
i=l 


satisfy  2 A  <  1.   We  shall  prove  in  section  5  that,  by  adding 
the  "stabilizing"  term 


K^T-h.1  V.  ut 
i=l 
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to  (4.1),  we  obtain  an  unconditionally  stable  difference  equation. 

It  is  not  difficult  to  verify  that  (4.2)  is  an  explicit  difference 
equation  when  the  net  points  in  a   are  suitably  ordered,  for  example, 
lexicographically . 

A  straightforward  application  of  Taylor's  theorem  shows  that  U 

satisfies  the  difference  equation  (4.2)  to  within  a  term  of  order 

9       N     , 
K+  |h|   +  K  \~~  h±x.      Therefore,  [11]  (4.2)  is  consistent  with  (4.1)  if 


i=l 

N    ,  -1 

k  =  o(i)  (y--  hi1)    , 

i=l 

as  n  — >  0. 

We  now  turn  to  the  alternating  direction  methods.   In  the  first 
alternating  direction  method  for  (4.1),  one  determines  an  admissible 
function  t  — >  u(t)   such  that 


V^tt)  -  u(t-K)]  =v1v1  u^)(t)  +yy1v1  u(t-K) 


K"1 


[u<1+1)(t)    -   u^(t)]    =Vi+1Vi+1   u(l+l)(t)    -V±V±   u(t-K), 


for      1=1,    2,...,    N-l,    and 


K-X[u(t)    -    u^-^Ct)]    =VN\j,u(t)    -     VNVNu(t-K), 

where  the  auxiliary  functions   t  ->  u*1' (t)  e  ^(T^),   defined  for 

t  e  S°,   are  such  that  u^1'  =  u   on  (^  X  S°.   We  assume,  of  course, 


-  13  - 


that   N  >  2,   otherwise  the  method  is  undefined. 

This  difference  scheme,  for  N  <  3,   was  investigated  in  [5], 
where  it  was  shown  that  the  operation  taking  u(t-K)   into  u(t) 
involves  the  inversion  of  p   tridiagonal  matrices,  where  p  is 
the  dimension  of   (^(fyj- 

Following  Douglas  and  Rachford,  we  eliminate  from  these  equations 
the  auxiliary  functions  u^  ' .      Adding  the  equations,  we  obtain 


(4.4)       u^-V^u  +  ^V^uW. 


From   the   last   equation 


(4.5)  u^"1)    =  u  -   K    VNVN  u  . 


Similarly,  solving  the  next  to  last  equation  for  u^    ' ,  we  find  that 
u(N"2)(t)  =  u(t)  -  K  Vn^Vn.!  u^^Ut)    +  K  V^V^  u(t-K), 

which,  in  view  of  (4.5),  becomes 

u(N"2)  =  u  -  K2  VNVN  u   +  Y?  V^Vn.^.!  u 


It  is  clear  that,  by  continuing  this  process  of  elimination,  we 
can  determine  the  functions   u^  '   as  a  linear  combination  of  u 
and  certain  of  its  difference  quotients.   When  the  resulting 
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expressions  for  the  u^  '      are  inserted  into   (4.4),  we  find  that 
u   is  a  solution  of  a  single  difference  equation  of  order  2N  +  1 


(4.6)        u_  =Ah(u)  +  y~   (-l)J+1KJ*DJ(u_), 
t  j=2  fc 


where  the  difference  operators   DJ   are  defined  by 


DJ  (u)  =  Y2*  Vi1Vi1Vi2  Vl2  * ' Vl  •  Vl 


J  "  "J 


and  the  sum  K  is  extended  over  the   (,)   different   j- tuples 

(i,  ,  i?,...,  i  .)   formed  from  the  first   N  positive  integers. 

The  usual  argument  involving  Taylor's  theorem  shows  that  U 

o 
satisfies  (4.6)  to  within  a  term  of  order  K  +  |h|  .   Hence  (4.6) 

is  consistent  with  (4.1). 

For  the  second  alternating  direction  method  we  restrict  ourselves 


to  the  case   N  =  2.   The  method  consists  in  determining  an  admissible 

,o 


function      t  — >  u(t)      such   that,    in     Jl     X  S, 


2K  2[v(t)    -    u(t-K)]    =V1V1    v(t)    +  V2V2   u(t-K) 
and 

2K_1[u(t)    -   v(t)]    =  V^   u(t)    +  V2V2   v(t), 

the  auxiliary  function   t  ->  v(t)  e  ^(i^),   defined  for   t  £  S°, 

»     o 
being  such  that  u  =  v   on  Jl  X  S„.. 
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Following  Peaceman  and  Rachford,  we  eliminate  the  function  v 
between  these  equations  and  find  that  u  satisfies  the  fourth 
order  difference  equation 

(4.7)       u_  =Ah(u)  -  f  Ah  (UJ  -  f  V^VgVg  u  . 

As  before,  we  see  that   U  satisfies  (4.7)  to  within  a  term 
of  order   K  +  |h|2. 

In  the  following  sections  we  shall  prove  that  these  difference 
equations  are  unconditionally  stable  in  the  following  sense:   if 
t  ->  u(t)  e  <§(Qh) ,      t  e  SR,   is  a  solution  of  (4.2),  (4.6),  or 
(4.7),  then 

l|u(t)||  <  C|u(0)||, 

N    -2 
where   C   depends  only  on  N   and  the  mesh  ratio   A  =  K(V~  h.  ). 

1=1 

5-   Stability  of  the  Semi-Explicit  Method. 

Before  stating  the  main  theorem  of  this  section,  we  prove 

several  lemmas.   These  will  be  seen  to  be  difference  analogues  of 

the  usual  quadratic  differential  identities  and  inequalities  which 

are  basic  to  the  energy  method. 

LEMMA  5 .   The  operators    V-   and   -  V-   are  adjoints  for  functions 

in   ^(ft^)?   that  is,  for  any  v   and   w  in   (?(£0 


(v,  V±  w)  =  -(Vj_  v,  w) 
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Proof.      We  have   the   identity 


(5-1)  vVi   w  =  V±[E  1(v)w]    -   E  ^-(^   v)    w. 


E|y    (2.1),      E"1(V1  v)    =S7±  v.      Since     v,    w  e   fc^) 


V.  [E_i(v)w]    =   0; 


and  we  conclude  from    (5.1)    that 


h 


N        r „„         ,.,     _  UN 


vyi   w  =  -   hIN     V~         (V±  v)   w, 


x   e    Gh  x  e    Qh 


which  is  equivalent  to  the  conclusion  of  the  lemma. 
LEMMA  4.   If   u  e  (^(fi^)   then 


(u,Ah(u))  -  -|Wli 


Proof.      Taking     v   =  u     and     w   =    SJ    u     in  lemma   3,    we  find   that 


(u,     V.V.    u)    =  -\\S7±  u||2, 


from  which  the  desired  result  follows,  by  summation  with  respect  to   i 
LEMMA  5.   For  any  function   t  ->  u(t)  e  ^(\) '  t   e  Sk'   —  have 
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2(u,    u.)    =    (||u||2)t   -   K[|uJ|2 


and 

2(u,    u_)    =    (||u||2)_  +  K   ||u_|p 
t  t  t 


Proof.   These  identities  are  immediate  consequences  of 


2  u  u   =  (u2)t  -  K(ut)2 


and 


2  u  u_  =  (u2)   +  K(uf_) 
t      T  t 

(See  [14]). 

LEMMA  6.   If   u  e  £({^)   then 

2(u,  ^h'1^  u)  „||u|g. 
1=1 

Proof.   From  lemma  J>,    with   v  =  w  =  u,   we  obtain  the  relation 

(u,  Vj_  u)  =  -(V±  u,  u). 

Since   V1  -  ^7^    =  h1  V.  V.*   we  conclude  from  this  that 
2(u,  V±   u)  =  -h±(u,  V±V±  u). 

Consequently, 
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2(u,  y-   h.  V±  u)  =  -(u,  Ah(u)) 

i=l 


=  IMg, 


by  lemma  4. 

LEMMA  7.   For  any  function   t  — >  u(t)  6  (^fl^)*  t  6  sR,   we  have 


2(ut,Ah  (u))  =  -(flu|g)t  -1-  K||ut||2. 


Proof.      From  lemma   3 


2(ut,V1Vi   u)    =  -2(V±   ut,v±   u) 


=  -  (l|V±   ul|2)t  +   K||Vi   ut||2, 


by  lemma  5-   The  result  follows  by  summation  over  i. 

N    p  1/2 
LEMMA  8.   Let   u.  =  f\~  h.  )    .   Then  for  any  function 

i=l 

t  -^  u(t)  £  (^(J^),  t  e  sR,   we  have 

2(u,   y--  h^1  V.    ut)    <   e    n  ||ut|r    +  i  H   ||u||J, 

i=l 


where   e  >  0  _is_  arbitrary, 
Proof.   From  lemma  3 


N    -1  N    -1 

(u,  y  h.   Vl  ut)  =  (ut,  y  h.    Vl  u) 


T=l  i=l 
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Applying  Schwarz'  inequality  to  the  right  side  of  this,  we  obtain 

(u,  y~  h"1  Vi  ut)  <  ||utn  \\y~  h"1  Vl  u|| 


i=i  i=i 


iV-llu^h-1!!^  u||  , 


i=i 


by  the  triangle  inequality  and  the  fact  that  ||V.  u||  =  ||V  u| 
Using  Schwarz'  inequality  again,  we  find  that 

N   -l- 
(u,  ^~  h.  V±  ufc)  <  ||ut||  M.||ul|x- 

1=1 
To  the  right  side  of  this  inequality  we  apply  the  generalized 
arithmetic-geometric  mean  inequality: 


(5.2)        2  ab  <  e  a2  +  |  b2       (e  >  o) , 


to  obtain  the  desired  inequality. 

THEOREM  1.   Let  the  function   t  -^   u(t)  e  ^,(^h),   t  e  S„,   be  a 

solution  of  the  semi-explicit  difference  equation  (4.2).   Then 

t  —  K 
l|u(t)||2  +  K^  ||u(n)||2  <  (1  +  2A(1  +  A))||u(0)||2, 

n=o 

which  implies  that  (4.2)  _is  unconditionally  stable. 

Proof.   The  difference  equation  (4.2)  is  satisfied  only  in  ft,  X  SK- 
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But  since  u(t)  e  £(£1  ),   u  vanishes  on  U^  X  SR,   and  hence 
we  may  form  the  inner  product  of  (4.2)  with  u  to  get 


(u,  ut)  =  (u,  Ah(u))  -  K(u,  y-  h.1  y.  ut), 


valid  for  t  e  S„.   From  this  and  lemmas  4,  5  and  8,  we  obtain, 

XV 

after  multiplication  by  2,    the   inequality 


(l|u|f)t   -    K||ut|f   +  2||u||2   <   e    K  M.||utir 
(5-5) 

If  we   choose      e      so   that      K  u.   =  e,      then    (5«3)    becomes 


(5-M  (||u||2)t    -    K(l   +   K  n2)||ut||2   +   ||u||2   <   0. 


Similarly,  we  may  form  the  inner  product  of  (4.2)  with 
u,   to  obtain 


In  veiw  of  lemmas  6  and  J,    this  becomes 


2iM2  -  -(HliV 
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When  this  result  is  combined  with  (5A),    we  obtain  the  inequality 


(IM|2)t  +|  (1  +  K  H2)(||u||2)t  +  ||u||2  <  0. 


Since 


t  —  K 
KY"   (l|u|p)t  =||u(t)||2  -||u(0 


rj=0 


we  have  from  this  that 


t-K 
|[u(t)||2  +  K  J"  ||u(T!)||2  <  ||u(0)||2  +  |  (1  +  A)||u(0)||2. 


r\=0 


o 
since  K  p.  =  A.   The  desired  result  follows  by  applying  lemma  2 

to  the  last  term. 

If  we  denote  by  e  the  error  U  -  u  for  the  semi-explicit 

method,  then  e(0)  =  0,   and  e(t)  e  £(%.)   for  t  e  SK-   Also, 

e   satisfies,  in  a  X  SK,   the  difference  equation 

N   -1„ 


et  -  Ah(e)  -  K^-h"  V±  efc  +  T, 

1=1 


where  the  function  t  — >  T(t)  e  ^(fyJ   is  the  local  truncation 
error.   It  satisfies 

t-K       0   1//2  0      N    -, 

(KY^IlTh)!!2)  =  0(K+  |h|2  +  K^T-h"1), 

T]=0  i=l 
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where  the  constant  implied  in  this  relation  is  independent  of 
K,   h   and   A. 

The  following  result  can  be  proved  along  the  lines  of 
theorem  1. 

THEOREM  2.   The  error   e  =  U  -  u  for  the  semi-explicit  method 
for  (4.1)  satisfies  the  inequality 

t-  —  K  t—  K 

||e(t)||2  +  K  JTH|e(T))||2  <  M  R  j-   |(  t(tj)||  2  ^ 

ri=0  ^=0 

2  P   N   -2 
where        M=^-+K+2lCy^  h,  . 
m  /I 


Proof.      As   in  the  proof   of   theorem  1,    we   obtain 

(5-5)        (I|e|!2)t-K|et|p  +  2||e|||  <  «   KuJkJ2  +  ±  K  |t||e|| 

+   2(e,T). 

In  view  of   Schwarz'    inequality  and    (5-2),    with     e      replaced 
by     £ ' ,    we  have 

2(e,T)    <   e'llefl2   +  A,    ||T|f. 

6 

According  to   lemma  1,      m|je||^   <  |le||^,      so   that 
2(e,T)   <|'||e|g  +  A,  ||t||2, 
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and  when  this  is  combined  with  (5«5)>  we  obtain 


K  e» 


e||t)*   -   K(l  +  €    Kn)||et|r   +    (2   -  ±  n  -  £    )||e||J 


<-  iiTir 


€» 


We  now  choose   e   and  e>   so  that  2  K  \i  =   £   and  2  e»  =  m, 

Then 


(5.6)  (||e||2)t   -   K(l   +  2A)||et||2  +||e|g   <§||T||2, 


As   before,    we  have 


l|et||2   =  -  |(||e||2)t  +    (et,    T) 


which,  after  invoking  Schwarz'  inequality,  becomes 


etH2  <  "  (He||2)t  +||  T||2, 


This  combined  with  (5-6)  gives 


e||2)t  +  K(l  +  2A)(||e||2)t   +||e||2   <    [§  +  K(l   +  2a)]||t||2 


from  which  the  desired  result  follows,  since  ||e(0)||  =||e(0)||-|   =  0. 

6.   Generalizations . 

We  now  indicate  briefly  how  the  results  of  the  preceding  section 
can  be  extended  to  parabolic  equations  of  the  form 
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(6.1)       i2.      y-_J_(a1J(x,t)  <*£-), 
it       ,  -C^-A-a,  A*, 


*'     i.fcr  <> 


for  x  e  ft  and  0  <  t  <  tQ.   We  assume  that  a1 ,  = 


a..   and  that, 


for  x  €  ft,   o  <  t  <  tQ, 


N 


1J 


p2ur  2    ;>_  *  J(*.t)e1eJ  >  pur, 


where   £  is  any  real   N-vector  and  p,  p,  >  0. 

We  approximate  (6.1)  by  the  semi- explicit  difference  equation 
(See  (2.2)) 


(6.2) 


N    -     i  i  N         lM  A 

ut  =  J2    V*   Vx  u)  +  y   ±(a1J  V.  u) 

i=l 


N  ,  _ 

-  K^   (p2  +  Pl  h"  )V±  ut, 

i=l 


where 


ii. 


a   (x,t)  =  a   (x,  , .  .  .  ,  x.  +  ?5 — , .  .  .  ,    x.j,  t) 


and 


p2  =    sup 
x  e  ft 

0  <  t  <  tQ 


ii 
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Denote  by  Lk(u)   the  sum  of  the  first  two  terms  on  the 
right  side  of  (6.2).   It  can  be  shown  that 

(6.5)        (u,  ^(u))  >  -  [P  +  0(|h|)].||u|j2, 

for  all  u  e  &(Oj . 

Similarly,  for  any  function   t  — >  u(t)  6  £(£1), 


(ut,  L^u))  <  (Au)t  +  P3l|u||^ 


(6.4) 


K 


f  (P2  +  |h|Pl)||ut|f  , 


where 


p3  = 


sup 
0  <  t  <  tQ 


da 

at 


ij 


and  _/\u  is  bounded,  both  above  and  below,  by  a  constant  multiple  of 

[1   +   0(K  +    |h|)]||u|g. 

With  the  aid  of  the  inequalities  (6. 3)  and  (6.4)  the  arguments 
employed  in  theorems  1  and  2  can  be  carried  over  to  prove  that  the 
semi-explicit  difference  equation  (6.2)  is  unconditionally  stable, 
provided  K  and   |h|   are  sufficiently  small. 
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7-   Stability  of  the  First  Alternating  Direction  Method. 
For   j  >  1  put 

1/2 

l|u||j  =  £^*ll?i  Vi  ..-Vi.  u||2) 

for  u  €  <^(a  ).   Then  we  have 
LEMMA  9-   If  u  e  £(\)      then 

N      1+1   1   1  N    1    ? 

(u,  y~  (-1)      K  D  (u))  -  -  T~  K  Hullr 

J^2  j=2 

Proof.   From  the  adjolntness  of  the  operators  ^7.   and   -y.   we 
have,  for  1  <  r  <  j, 


(u,  Vi  vln  •  •  -Vi  .Vi  .  u)  =  "(Vi    u>  Vi  vin  •  •  -Vi  •  •  -Vi  .  Vi  .  u) 

11  JJ  r  11  r  JJ 


From  this  we  obtain  the  relation 


1   1    j   J*  1     j 

Summing  both   sides   of   this   relation  over   the      (,)      j-tuples,    we 
find  that 


(u,    DJ'(u))    =    (-l)J  ||u||2. 


The  desired  result  follows  from  this  by  multiplication  by 
(-1)J+   KJ   and  summation  over   j. 
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The  arguments  of  lemmas  7  and  10  can  be  employed  to  prove  the 
LEMMA  10.   For  any  function  t  — >  u(t)  6  £(fyu) ,  t  6  SR,   we  have 

2(u,  y~    (-1)J'+1  KJ  DJu_)  =  -  y~   KJ(||u||2)_ 

3=2  t       3=2  t 

-y~   KJ'+1||u_||2. 
3=2  t 

THEOREM  5-   If  the  function  t  ->  u(t)  6  C*^) ,   t  e  sR,   is  a 
solution  of  the  difference  equation  (4.6) ,  then 

||u(t)||2  +  2K^||u(t1)||2   <  c||u(0)||2, 

T)=K 


where 


c  =  (1  +  4a)N  -  4*. 


Hence,  the  first  alternating  direction  method  is  unconditionally  stable. 
Proof.   Since  u(t)  e  £(^0  for  each  t  €  s°,   we  have 

(u,  u  )  =  (u,Ah  u)  +  (u,  y~   (-l)j+1  Kj  Dj  u  ), 

fe  j=2  fe 

valid  for  t  e  s£.   Using  lemmas  4,  5,  and  10,  this  becomes 

is. 

(||u||2)_  +   K||u_||2   +   2|u||f   +y~  KJ'(||u||2)_ 

t  t  fa  J  t 

=  -  J"  K^  ||u    ||2. 
j=2  t 
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Dropping  the  second  terra  on  the  left  and  the  negative  term  on  the 
right ,  we  find  that 

(7.2)      ||u(t)||2  +  2K  Y~   ||u(n)||2  <  ||u(0)||2  +  Y2   KJ  Hu(°)ll5- 

r|=K  j=2 

But,  from  lemma  2,  we  find  that 


KJ  ||u(0)||3  =  KJ  y~  ||V±  ...V±   u(0 


)ll 

J 


J  ir J  11,,/nMI2  ^ *  k"2    h~2 


<  4J  KJ  II U 


cor  y~  h"-=...h-' 

^—      i 


<  4J  l|u(0)rAJ  (?), 


where 


A  =  max  Kh~j  , 
i    J 


Therefore, 


KJ'  ||u(0)||2  <  [(1  +  4A)N  -  4nA-  1]  ||u(0)||2 


j=2 


[(1  +  4a)n  -  1  -  4a]  ||u(0)||2, 


and  this,  together  with  (7.2),  completes  the  proof  of  the  theorem, 

As  before,  we  can  prove  the 
THEOREM  4.   If   e  =  U  -  u   is  the  error  for  the  first  alternating 
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direction  method,  then 


l|e(t)||2  +  K^||e(n)||2  <IK^||T(T,; 

Tj=K  ti=K 


where   t  — >  T(t)  £  (^(ty-)   is  the  local  truncation  error,  and 
T  -  0(K  +  |h|2). 

Similar  results  can  be  proved  for  parabolic  equations  of 
the  form 


(7.3)       <^U  =  y~  ±-     (alj'(x,t)  ^  ). 


3t   frj^i  <)x± 

8.    Stability  of  the  Second  Alternating  Direction  Method. 

LEMMA  11.   For  any  function  t  — >  u(t)  e  (^(o^) ,  t  e  SR,   we  have 

2(u_,Ah  u)  =  -  (||u||2)_  -  K||uJ|f 

t  t       t 

The  proof  is  similar  to  that  of  lemma  7  and  will  therefore  be  omitted. 
THEOREM  5-   If  the  function   t  ->  u(t)  e  <?(^) ,   t  e  SK,   is  a 
solution  of  the  difference  equation  (4„7),  then 

||u(t)||2  +  2kyH|u(t1)||2  <   [H  -  8a  +   (1  +  4a)2]  ||u(0)||2, 


where   C  =  A  if   4a  <  1   and  C,   =   4a  -  1/2      if   4a  >  1.   Hence, 
the  second  alternating  direction  method  is  unconditionally  stable, 
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Proof.   From  (4.7)  we  have,  as  before,  that 


(u,  u_)  =  (u,  \   u)  -  I  (u,  ^   u_) 

t  w 


K2     -   - 

-  ^  (u,  V1V1V2V2  u_) 


In  view  of  lemmas  4,  5,  and  10,  this  becomes 


u||2)_  +   KlIuJI2  +  2[|u|g   =|   (||u||2)_ 
t  t  t 


+  f2  II  u    l||   -^    (||u|||)      -$5||u_|||. 


Taking  the  inner  product  of  (4.7)  with  u_,   we  find  that 

t 

(8.2)        ||u_i|2  =  -  |  (||u||2)_  -  f^ujg. 
t  t       t 


We  now  multiply  (8.2)  through  by  2  C,   K  and  add  the  result 
to  (8.1)  to  obtain 

(||u||2)      +  K(2    C  +   DlluJI2   +  2||u||2  +    (C  -  |  K(||u||2)_ 
t  t  t 

(8.3) 

<|Wlli  -  -f  (IM|)_. 
^      t  -1-  t 
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According  to  lemma  2, 


K^ujg  <  4Aiiu_u2, 

t  x  t 


so  that  (8.3)  implies  that 


u||2)      +    (2C  +  1   -    ^)K||u_||2   +   2||u||2   +    (C   -  |)K(||u||2)_ 

t  t         x        *      a  t 

+  f(\W\\22)_<  0. 


In  view  of  the  definition  of   £,   we  have  from  this  that 
||u(t)||2  +  2K^|[u(t1)[|2  <  ||u(0)||2  +  (C  -  |)K||u(0)«2 

T)=K 

+  fWo)||2, 

and  this,  with  the  aid  of  lemma  2,  gives  the  desired  result. 

Similarly,  we  have  the 
THEOREM  6.   If   e  =  U-u  _i£  the  error  for  the  second  alternating 
direction  method,  then 

||e(t)||2  +  K^||e(n)||2  <  (|+  C)K^||T(n)||2, 

T]=K  T]=K 

where   £  _is  defined  in  theorem  5,  and   t  ->  T(t)  e  &(%)   is  the 

o 
local  truncation  error;   T  =  0(K  +  |h|  ) . 


Similar  results  can  be  proved  for  the  parabolic  equation  (7-3) 
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